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1.  The  Pioblem  , v 

; \ 

A my  may  he  defined  as  the  locus  of  piopagatlon  of  an  element  of 
wavefront.  The  definition  includes  the  plane-wave  idealization  that  the 
wavefront  propagates  along  Its  Instantaneous  normal;  a propagating 
wavefront  is  discretized  by  a selection  of  lays;  each  lay  Is  bent  by 
infraction  in  a medium  having  sound-speed  gradient.  My  definition,  lays 
me  oithogonal  to  wavefronts. 

We  denote  the  specific  wavefront  length,  in  a vertical  plane,  by 
l,.  Along  a lay  as  the  wuvetiont  diverges  l.  Increases;  as  the  wavefront 
conveiges  1.  decreases  and  may  pass  through  zero.  Passing  through  zeio 
may  Indicate  a folding  oi  focusing  ot  the  wavefront.  A eaustU  Is  a ray-to- 
ray  locus  ot  an  l.  * 0 occurance.  A caustic  generally  indicates  a folding  ^ 

ot  the  wavetront.  The  end  points  ot  a caustic  may  be  indicative  of 
focusing . 

We  specialize  to  the  case  of  sound-speed  varying  with  depth,  z, 
only,  and  a horizontal  caustic  at  depth  zA.  This  Imposes  horizontal 

uniformity,  without  caustic  end  points.  \ 

l oi  the  analysis  we  assume  tin'  sound-speed  profile 

c - U z (1) 

where  p Is  the  gradient  at  caustic  depth,  and  the  zero  depth  level  is 

defined  by  the  speed,  c^ , at  the  caustic  depth.  We  shall  have  to  decide 

latei  as  to  the  extent  of  the  region  of  analysis  applicability  about  the 

♦ 

caustic. 

The  associated,  consistent,  ray,  wavefront  and  caustic 
configuration  Is  depicted  in  Fig.  1. 
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According  to  ray  theory  the  spreading  factor,  amplitude  A,  for  sound 
intensity  is  given  by 


A 


(2) 


Sign  changes  in  L along  a ray  are  indicative  of  wavefront  foldings. 

At  a caustic  L = 0 and  A -»  oo  . Does  this  mean  that  every  point 
on  a caustic  is  a focal  point?  Generally,  effectively  not.  We  have  both 
incoming  rays  and  their  refracted  return  (except  at  caustic  ends)  and  it  is 
the  resultant  that  manifests.  We  here  attempt  to  determine  the  resultant 
by  analysis  consistent  with  ray  theory.  That  is,  we  wish  to  extend  ray 
tracing  for  applicability  at  caustics. 
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2.  The  Analysis 


For  a linear  speed  profile  the  rays  are  all  circle  arcs  with  centers 
at  z = 0 and  radius  z^  * • ^or  these  rays 


cos  V = — -= 


sin  y ■ + 


r 2 2 
£ z*  " 2 1 


We  here  adopt  the  convention  that,  should  there  be  a sign  duality,  the 
upper  sign  applies  to  incoming  rays/wavefronts  and  the  lower  sign  to 
returning  rays/wavefronts.  The  rays  may  be  expressed  by 


- z#  sin  y 


where  x is  introduced  as  the  horizontal  linear  coordinate.  Each  ray  is 
identified  by  x^;  the  circle  arc  center  is  at  x#  and  z = 0,  and  the  ray 
touches  the  caustic  at  x^,  z^  . 

The  specific  wavefront  length  may  be  expressed  by 


6x^  sin  y 


And,  by  ray  theory,  the  normalized  wave  amplitude  may  be  expressed  by 


A * ± nrr 


3- 


The  wavefronts  are  obtained  by  integration  of 

T*  tan?  (») 

dz 

which  deilves  from  theti  being  orthogonal  to  the  rays.  The  wavefronts  may 
be  exptessed  by 


**  sin  ? 7 ** 


cos  y 


X - X , 


<!») 


Kach  wavefront  (line  of  constant  phase)  is  Identified  by  x#l  where  the 
wavetiont  touches  the  caustic. 

At  tin1  caustic  limit  the  wavefronts  are  moving  hot izontally  with 
sped!  c#  . A steady  state  conligui  ation  is  realized  front  a trame  of 
toferenoe  moving  with  that  velocity.  Wi'  ailopt  this  frame  of  reference. 

At  this  point  we  Introduce  the  waveform 

C>  ■ A cos  kxA  (10) 


which  corresponds  to  the  lieguency 

U)  k c#  / 2tr 


(U) 


Hy  substituting  f oi  x^  from  Kg.  (9)  and  for  A from  Kg.  (7)  we 
may  exptess  the  incoming  wave  by 

O.  - ~“r  cos  k [x  - B (z)l  (12) 

i sin  y i 
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where 


e.(z) 


- z^  sin  y + z*  Bm 


1 + sin  y 
cos  y 


(13) 


And  the  returning  wave  by 


O 

r 


1 

sin  y 


cos  k [x  - (z)] 


(14) 


where 


(z) 


- z^  sin  y 


z^  Sn 

★ 


1 - sin  y 

cos  y 


(15) 


For  the  resultant,  4>  , we  add  the  two  fields: 

$ = o.  + O (16) 

i r 

In  expressing  the  resultant  we  use  positive  y and  sin  y only.  Thus 
for  the  returning  wave  the  signs  of  sin  y are  reversed  in  performing  the 
addition: 


* 


cos  k [x-0]  - . 1 'C  cos  k [x+0] 


sin  y 


sin  y 


sin  y 


sin  k 0 sin  k x 


(17) 


where 


’tn  i-± 
c< 


sin  y 


cos  y 


- sin  y 


t 


(18) 
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Lot  us  examine  the  resultant  in  tiro  region  of  the  caustic,  whore 


y <<  1.  Expansion  reveals  that  0 may  be  approximated  by 

6 * z*  j V3  (19) 

For  the  supplemental  condition  that  k 0 < <•  1,  the  resultant  may  be 
approximated  by 


•I  tt  re  * 

G'  ■*>  — ~ y ~ sin  kx  (20) 

where  we  have  replaced  kz*  according  to  Eqs.  (1)  and  (11).  From 
Eq.  (-1)  we  derive  for  small  V,  small  6 z z^  - z , that 

> " « 6z  (21) 

c* 

Equations  (20)  and  (21)  combiire  to  reveal  that  the  resultant  approaches  a 
flrst-ordei  zero  at  the  caustic: 


<I> 


OJT 

3 


<*L 

C. 


6 z sirr  kx 


(22) 


By  ray  theory  the  spreading  factor,  amplitude  A,  is  zero  at  a caustic. 
This  analysis,  however,  does  irot  apply  in  the  region  of  caustic  end 
points . 

It  is  relevant  to  determine  the  depth  and  value  of  maximum 
amplitude  for  the  resultant,  G1 , provided  we  can  show  that  it  occurs  near 
the  caustic  in  the  region  of  analysts  applicability.  We  denote  the  level 


by  z and  values  there  by  subscript  m.  We  first  assume  that  y<<  1 
m 

and  later  show  this  to  be  so. 

For  small  V we  approximate  0 by  Eq.  (19)  and  express  the 
amplitude  factor  of  Eq.  (17)  by 


A = 


2_ 

y 


w 


y 


3 


The  extremal  occurs  where  dA/dV  = 0.  This  yields 


tan 

2n  u 3 

= 3 

i 

CO 

31 

£|< 

3pm 

p m j 

M w y 3 
3 p m 


* 1.32 


y * 

1.32  y- 

U 

1/3 

» 0.86 

m 

2n 

U) 

U! 

1/3 


A 

m 


y 

m 


2.26 


1/3 


(2  3) 


(24) 


(25) 


(26) 


(27) 


In  order  to  evaluate  V we  must  specify  p/u)  in  Eq.  (26).  Even  for  the 

m ^ 

relatively  high  value  of  10  we  obtain  y m *>  0.086  which  satisfies  our 

condition  of  smallness;  corresponding  thereto  A ^ 22.6,  a considerable 

m 

amplification. 

Equation  (25)  expresses  the  phase  angle  k©m  which  appears  in 

O and  0 as  combined  for  Eq.  (17).  We  note  that  z occurs  somewhat 
i r m 

below  k 0 = v/l  which  is  the  level  where  <p  , 0^,  and  4>  are  all  in 
phase. 
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From  Eq.  (3)  we  obtain 


* 

z = — cos  y (28) 

m /i  m 

For  c = 1500  m sec  n = 10  1 sec  \ and  ui  = 100  Hz  giving 

y = 0.086,  we  obtain 
m 

z.  - z 83  54  m (29) 

* m 

Analysis  applicability  requires  that  the  speed  profile  be  linear  about  the 
caustic.  This  linearity  must  extend  to  beyond  the  z^  level  for  the 
maximum-amplitude  analysis  to  be  valid. 


3.  The  Result 

The  payoff  of  the  analysis  lies  in  Eq.  (27).  We  can  split  the 
resultant  amplitude , <I>,  between  o.  and  <p  to  get  at  an  effective  lower 
bound  on  L in  computing  the  amplitude  factor  along  a ray.  The  analysis 
suggests  the  approximate  lower  bound: 


L * 
m 


(30) 


However  we  must  not  forget  the  normalizing  factor.  Equation  (7)  indicates 
how  we  have  normalized  L.  We  now  write 


L = sin  (y  - yQ) 


(31) 
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where  7 is  the  ray  angle  at  caustic  tangency;  that  is,  we  do  not  require 
the  caustic  to  be  strictly  horizontal.  The  missing  factor  in  L of  Eq.  (31), 
lost  by  normalizing,  can  be  obtained  from  L , along  the  ray,  before  and/or 
after  the  ray  passes  through  the  caustic.  Adding  this  factor  to  Eq.  (30) 
yields 


L 

m 


(32) 


t 


This  is  where  ray  theory  has  led  us.  The  approximation  omits 
caustic  curvature  and  lateral  variability  in  rays.  Equation  (32)  also  does 
not  apply  near  caustic  end  (branch)  points  where  interference  is  partial  or 
nil  and  strong  focusing  may  occur.  It  remains  to  reconcile  the  results 
with  full  wave  theory  and  energy  considerations. 
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